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ABSTRACT 

In a previous paper (Falkenberg, KotuUa & Fritze 2009, astro-ph/0901.1665) we have 
shown that the classical definition of E+A galaxies excludes a significant number 
of post-starburst galaxies. We suggested that analysing broad-band spectral energy 
distributions (SEDs) is a more comprehensive method to select and distinguish post- 
starburst galaxies than the classical definition of measuring equivalent widths of (US) 
and [O ii] lines. 

In this paper we will carefully investigate this new method and evaluate it by 
comparing our model grid of post-starburst galaxies to observed E+A galaxies from 
the MORPHS catalog. 

In a first part we investigate the UV/U-optical-NIR SEDs of a large variety in 
terms of progenitor galaxies, burst strengths and timescales of post-starburst models 
and compare them to undisturbed spiral, SO and E galaxies as well as to galaxies in 
their starburst phase. In a second part we compare our post-starburst models with 
the observed E+A galaxies in terms of Lick indices, luminosities and colours. We then 
use the new method of comparing the model SEDs with SEDs of the observed E+A 
galaxies. 

We find that the post-starburst models can be distinguished from undisturbed 
spiral, SO and E galaxies and galaxies in their starburst phase on the basis of their 
SEDs. It is even possible to distinguish most of the different post-starburst by their 
SEDs. From the comparison with observations we find that all observed E+A galaxies 
from the MORPHS catalog can be matched by our models. However only models 
with short decline timescales for the star formation rate are possible scenarios for 
the observed E+A galaxies in agreement with our results from the first paper (see 
Falkenberg, KotuUa & Fritze 2009a). 

Key words: Galaxies: evolution - Galaxies: formation - Galaxies: interactions - 
Galaxies: starburst - Galaxies: clusters: general 



1 INTRODUCTION 



The morphology-density relation (lOemleij Il974l: iDressler 
1980*) and the Bu tchor-Ocmler effect ( Butcher fc Oemler 
'l978, 19ii;| Yang et al. 20041') indicate a transformation from 
late- to early-type galaxies during the lifetime of a galaxy 
in high density environments. Over the past a wide range of 
transformation processes have been proposed, acting in dif- 
ferent environmental densities and on different timescales: 
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Important mechanisms are mergers, harassment, gas strip- 
ping and strangulation. 

Merging of two or more galaxies is an interaction most 
efficient at low relative velocities. If at least one gas-rich 
galaxy is involved in the merger this triggers a starburst 
and a subsequent morphologica l transformation, in genera l 
into an early- type g alaxy (e.g. Barnes fc Hernauis1ill996l ). 
Galaxy harassment (|Moore et al.lll996l . f 19981 ) are multiple 
high-speed encounters in the densest regions o f the cluster, 
but c an have effects out to large cluster radii (|Moran et al.l 
l2007f ). This harassment proc ess can also be connected to 
a starburst (|Lake et al.lll99^ ). Finally there are gas strip- 
ping and strangulation, two related processes in which the 
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infall into the galcixy cluster and contact to the hot X-ray 
gas within the cluster leads to a truncation or termination 
of star formation (SF). One example is ram pressure strip- 
ping (RPS) of the infalling galaxy's gas. The difference be- 
tween truncation and termination is the timescale on which 
the gas is removed: In the case of truncation the RPS is 
strong enough to strip the gas out of the disk, truncating 
SF on timescales of ~ 100 Myr. Termination on the other 
side primarily strips off the halo gas, leading to a halting 
of SF on longer timescales of ~ 1 Gyr dKawata fc Mulchaev 



20011 ) 



20081 : iMoran et aUbOOd : iBarr et al.ll2007l : lKodama fc Bower 



To evaluate the possible mechanisms that take place 
during this transformation in cluster environments like ha- 
rassment, gasstripping, strangulation and mergers (see, e.g., 
|Poggiantill2004al ). the investigation of post-starburst galaxies 
as intermediate state between late- and early-type galaxies 
is necessary. Several ea rlier studies (e.g.. iBarger et al]|l996l : 
iLeonardi fc Rose! 1 19961 ) found that a significant fraction of 
cluster members show spectral features characteristic for a 
recent starburst in the last 2 Gyr, and that those starbursts 
affected large fractions of the total galaxy mass. If those star- 
bursts are triggered by merging this can explain the transfor- 
mation of disk-dominated into bulge-dominated (SO) galax- 
ies. This would naturally explain the resulting decrease in 
the fraction of blue galaxies in interme diate rcdshif t clusters 
as compared to local sampl es (|Dressler,.1980. : ,Drcssler et al.l 
1 19971 : IPoggiaiit i et al.ll2008l ). An alternative explanation for 
the same effect is the truncation of SF in spirals and the 
subsequent passive evolutio n, leading to a fading and migra- 
tion o nto the red sequence (|Treu et al.ll2003l : [Pe Lucia et al.l 
|2004 ). 

iDressler et"aL I (|l999l . Table 6) introduced a classifica- 
tion scheme based on the equivalent widths of [O ll] and H5 
lines. The relevant classes for the present paper are k, k-l-a, 
a-|-k and e(a). e(a) galaxies have spectra that have strong 
Balmer absorption EW(H(5) > 8 A but [O ll] in emission. 
The remaining classes a-l-k, k-f a, and k all have no deteca- 
ble [O ll] emission, but differ in the strength of their balmer 
absorption, with a-l-k having strong (EW(II(5) ^ 8 A), k-|-a 
moderate (sA ^ EW(H5) sA) and k-type spectra weak 
(EW(II(5) ^ 3AA) KW features. Another common term for 
post-starburst galaxies is the class of E-l-A galaxies. E-l-A 
galaxies are described by a superposition of an elliptical 
galaxy spectrum (hence the E) with the spectra of A-type 
stars with strong Balmer absorption lines (the A in E+A), 
making if a more general class including a-l-k and k-|-a 
galaxies. In the remainder of this paper we will therefore 
use the term E-l-A interchangeably with the combination 
of k-|-a/a-fk galaxies. Note that we assign positive signs to 
EW(H5) if H(5 is in absorption, while positive signs mean 
emission of the [O ll] lines. 

In a first paper jPalkenberg fc Fritzel |2009| . IPaper 1 
hereafter) we focused on modelling E-l-A galaxies (|Gotc 
l2004l : IPoggianti et al.|[l999l ) with our galaxy evolution code 
GALEV. For this purpose we imposed different combina- 
tions of a starburst and halting of star formation (SF) with 
different timescales and beginning times on spiral galaxy 
models. 

We found that the classical definition of E-l-A galaxies 
by measuring the eq uivalent width of Balmer absorption and 
[O ll] emission lines (|Gotdl2004l : |Poggiantill2004bl : lYang et all 



|2004 ) excludes a significant number of post-starburst galax- 
ies, in particular those in stages with still some SF going 
on after a burst. By investigating model spectra we predict 
that Spectral Energy Distributions (SEDs) extending from 
UV/U through optical and eventually near- infrared (NIR) 
passbands 1.) allow a better distinction of E-|-A progeni- 
tors and scenarios and 2.) give a comprehensive census of 
all post-starburst galaxies of which, as we have shown, the 
E+As are only a subclass. 

In the first part of this paper we investigate the full 
SEDs of our models. The second part shows a comparison 
of models with observations from the MORPHS group. 



2 GALEV MODELS 

To model the galaxy transformation pr ocesses we use our 
GALEV evolutionary synthesis models ( Bicker et al.l |2004| : 
lAnders fc Fritzdl2003l : iKotuUa et al.ll2009i ). The full descrip- 
tion of the used models and the modifications to compute 
the [O ll] emiss ion line strength and H(5 equivalent widths 
can be found in |Paper~3 . hence we only give a short overview 
here. GALEV models allow to trace the spectral and photo- 
metric evolution of galaxies with arbitrary star formation 
histories (SFH s). We use stellar evolution data from the 
Padova group (iBertelh et al.ll 19941 ) including the important 
TP-AGB phase. To accurately model gaseous emission lines 
we supplement these d ata with data for the Z ero Age Main 
Sequence (ZAMS, see lAnders fc Fritzel l2003l, for details) 



We fu r ther u se the stellar template library of 



Leieune et al 



1997, 1 19981 ). Convolving the isochrones with a Salpeter 



(1955D-IMF with mass-limits of 0.1 and 100 Mq and steUar 
template spectra yields isochrone spectra that can then be 
weighted with the SFH of each galaxy type and integrated to 
yield the galaxy spectrum. Magnitudes in the various bands 
are then computed by convolving the galaxy spectra with 
the filter response curves and applying the appropriate zero 
points. 



2.1 Undisturbed galaxies 

SFHs of undisturbed spectral types are calibrated using 
a wealth of observational data such as integrated spec- 
tra and colours, star formation rates (SFRs), gas-fractions, 
gase ous and stel l ar me tallic ities, and mass-to-li ght ratios 
(see iBicker et all (|2004l ) and iKotuUa et al] (|2009l ) for a de- 
tailed description of this calibration). Note that our galaxy 
types E/SO and Sa to Sd are meant to represent spectral 
types resulting from the chosen SFH. In the local universe 
those spectral types correlate very well with morphological 
types, but this might not hold to arbitrarily high redshifts. 
To fully reproduce the observed properties of local, observed 
galaxies it is furthermore necessary to match the observed 
absolute magnitudes for the average representative of each 
galaxy type. Since colours and spectra only depend on the 
SFH (which is assumed to be independent of galaxy mass) 
we chose a total galaxy mass (including both stars and gas) 
that matches after a Hubble time the average absolute mag- 
nitudes of galaxies from the Virgo galaxy cluster. 
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Table 1. Parameter space covered by our grid. The upper part 
lists the models we use to demonstrate the effects of different 
starbursts scenarios on the observed SEDs, the lower part gives 
details on the models we compare to observations (Scct.|4]l. The 
latter also include evolutionary and cosmological corrections. 
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2.2 Modelling galaxy transformation scenarios 

To model galaxy transformation scenarios we add a star- 
burst followed be a truncation of SFR or a pure trunca- 
tion of SF dependi ng on the transformation process (see 
Sect. 2.4 of IPaper 11 for details). Starbursts are described by 
three parameters: 1) the time tburst or redshift Zburst of the 
burst; 2) the burst strength 6, i.e. the fraction of gas avail- 
able at the onset of the burst that is converted into stars 
during the burst; 3) the decline timescale r of the exponen- 
tial declining SFR after the burst. The full parameter set 
used in the rest of this paper are summarized in Table [T] 
Note that this model grid does not cover the full parameter 
grid. The reason for that is that we do not list any mod- 
els that do not fulfill the E-l-A criterion (e.g. models with 
Sa progenitors and/or weak bursts). In this work we iden- 
tify galaxies as k-l-a type galaxies if they have H5 equivalent 
widths EW(H5) ^ 5 A and only weak [O ll] emission lines 
with EW([Oll]) < 5A. 

2.3 Impact of dust extinction 

Besides the SFHs the amount of dust extinction is the 
dominating factor influencing the colours of galaxies. How- 
ever, in the present study we are not dealing with galax- 
ies with ongoing starbursts, but rather galaxies in their 
later, post -starburs t phases with no or only little remain- 
ing SF. In IPapei^ we found that only galaxies with strong 
bursts 6 ~ 70% reach the H5 strong phase. By the time 
it reaches this stage the galaxy should be essentially gas- 
free, since most of the gas was consumed in the past burst 
and the remaining gas ejected due to feedback processes 
(see Sect.[5ll. This is also confirmed by non-det ection of the 
majority of k-l-a gal axies in HI observations l|Gotd |2004| : 
iMiller fc OwenI l200ll 'l. Since the amount of dust is linked 
to the gas reservoir we do not expect significant amounts 
of dust in t hose late post - starbu rst phases. This is also 
supported by iBalogh et al.l (|2005h who used dust sensitive 
optical-NIR colours and found that their spectra are not sig- 
nificantly affected by dust. Note that this does not contra- 



Table 2. 1 a ranges for o b served galaxies in the Virgo cluster, 
taken from ISandaee et al.l lll985bl . Fig. 10) and ISandage et al.l 
lll985al . Fig. 4) 
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diet the findings of (very) dusty ongoing starbursts during 
earlier phases of ga l axy transformation , e.g. the e(a) phases 
JShioya ct al."200l'; iBalogh et aLll2005l : jPressler et al I I2OO8I : 
IPoggianti ct al..,200i '). 



3 SPECTRAL ENERGY DISTRIBUTIONS OF 
POST-STARBURST GALAXIES 

With SEDs it is possible to investigate the changes in mag- 
nitude for each filter as well as in colours of a galaxy at 
the same time. This enables an easy comparison of SEDs 
of undisturbed galaxies with SEDs of galaxies undergoing a 
starburst or SF truncation event. 

In Fig. [T] we show the SEDs of undisturbed Sa and Sd 
models. The symbols show magnitudes in t he Vega system 
for t he Johnson U to I and Bessell-Brett (jBessell fc Brett! 
Il988l ) J to K filters, put at their respective central wave- 
lengths. The lines connecting the symbols are only meant to 
guide the eye. We recall that the luminosities of our undis- 
turbed model galaxies at an age of 13 Gyr are gauged to 
agree with average observed B-band luminosities of the re- 
spective galaxy types Sa, Sb, Sc and Sd as observed in the 
Virgo cluster. However, real galaxies show a spread in ab- 
solute magnitudes. The observed 1 a ranges for galaxies of 
each type, derive d from luminosity func tions of galaxies in 
the Virgo cluster (|Sandage et al.|[l985bl lah. are given in Ta- 
ble [3 

3.1 Luminosity evolution 

Galaxies become fainter and bluer towards later types. Mag- 
nitudes of the typical Sa-type galaxy (see top panel of Fig.[l]) 
hence range from ~ —19 mag in the U-band to ~ —23 mag 
in the K-band. The typical Sd-galaxy in comparison (bot- 
tom panel of Fig. [1} is fainter by ~ 2 mag in the U-band and 
~ 3.5 mag in K. 

The time evolution of typical brightnesses differs be- 
tween galaxy types. Early-type spirals, such as Sa-type 
galaxies, with their declining star formation rates (SFRs) 
are bright early in their life and become fainter with time. 
This effect is most dramatic in the U-band because short 
wavelengths are most affected by the lower SFR at later 
times, while long wavelengths show smaller changes. For Sd 
galaxies and their constant SFR scenario the blue part of 
the spectrum, which is dominated by massive, hot stars with 
short lifetimes, very quickly reaches an equilibrium state and 
remains constant with time. However, cooler, low-mass and 
hence long-living stars slowly accumulate with time, leading 



4 Falkenberg, Kotulla & Fritze 



-23 



-22 



F -21 



^ -20 



-19 



-20 



-19 



— •— 5.0 Gyr 

- — 6.0 Gyr 
■ 6.5 Gyr 
□ 7.0 Gyr 
8.0 Gyr 
- -o- ■ 10.0 Gyr 






- 


































/uYb \ /vYr 








4000 6000 


8000 10000 
Wavelength [A] 


16000 


22000 



■■^ -IE 



^ -17 



-16 



5.0 Gyr 
6.0 Gyr 
6.5 Gyr 
7.0 Gyr 
8.0 Gyr 
10.0 Gyr 




4000 6000 8000 10000 16000 

Wavelength [A] 



22000 



Figure 1. SEDs for an undisturbed Sa (top panel) and Sd model 
(bottom panel) with ages between 5 and 10 Gyr. The connecting 
lines are only meant to guide the eye. 



to an increase of luminosity with time in the long- wavelength 
bands. 

Undisturbed Sb and Sc models evolve similar to the Sa 
model, while the luminosity difference between the different 
ages is smaller, i.e. the SEDs for different ages lie closer 
together. Luminosities of the Sb and Sc models vary between 
those of the Sa and Sd models. 

Note the change with time in the shape of the SED at 
the short wavelength end in case of the Sa model with its 
strongly decreasing star formation rate as opposed to the 
case of the constant SFR Sd model. 

Fig. [2] shows the SEDs of an Sd model undergoing a 
burst leading to the maximum peak value for EW(H^ ) and 
the longest E-l-A phase of our model grid (see lPaper J ). The 
burst starts at a galaxy age of tburst = 6 Gyr and is de- 
scribed by a sudden increase of the SFR to a peak value and 
a subsequent exponential decline of SFR with time on a char- 
acteristic timescale r, here r = 10* yr. The burst strength, 
in this case 70%, is defined as the fraction of gas available 
at the onset of the burst and converted into stars during the 
burst. 

At 5 Gyr the model is still undisturbed, the SED identi- 
cal with the SED at 5 Gyr in the bottom panel of Fig.[T] At 6 
Gyr all luminosities of the burst SED have increased and the 
galaxy is now much brighter than the undisturbed model at 
any time. The overall shape of the SED has changed, i.e. the 
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Figure 2. SED of an Sd burst model with ages between 5 and 
10 Gyr (Sd, burst of 70% at a galaxy age of 6 Gyr and decline 
time T = 0.1 Gyr). 



starburst galaxy model is even bluer than the undisturbed 
Sd-model. U-B colours, which appear as kink between U and 
B in SEDs like the one in Fig. (2] can reach —1 mag for the 
strongest bursts. This kink is, although less prominent, also 
present in SEDs of undisturbed spiral galaxies, but absent 
in SEDs of nearby SO and E galaxies. 

At an age of 6.5 Gyr, 0.5 Gyrs after the burst, the 
galaxy already is in its H5-strong phase, but has still mea- 
surable [O ll] lines. We recall that we consider a threshold 
of EW([Oll])= SA; only galaxies with a lower equivalent 
width are regarded as E-l-A galaxy with absent [O ll] lines. 
This threshold of 5 A is t he upper t hreshold for EW([Oll]) 
from iDressler et al.l (|l999l ) (see also lPaper j ). Since the be- 
ginning of the burst luminosities in I to K have increased, 
while the U to R luminosities already started to decrease, 
but are still brighter than the undisturbed model before the 
burst. The H5-strong model is now redder than the starburst 
and undisturbed Sd models, the SED of the undisturbed Sd 
model. 

The SED 1 Gyr after the burst shows the typical SED 
of an E-l-A galaxy, for which the EW([Oll]) is now lower 
than 5 A. For the SEDs at 7, 8, and 10 Gyr, all luminosities 
decrease with increasing age. The shape of the SED evolves 
with time into the SED of an early-type galaxy, which is 
even redder than a galaxy in the H5-strong and E+A phase. 



3.2 Impact of burst parameters on the luminosity 
evolution 

In general the increase in luminosities and change in colours 
during the burst and the entire post-starburst phase, i.e. the 
phase in which the SEDs can be distinguished from those of 
the undisturbed spiral galaxy on the one hand, and from the 
SO/E galaxy SEDs on the other, is influenced by the impact 
of the burst. This impact in turn depends on burst strength 
b, decline time r, onset of the burst, and galaxy type. 

Progenitor galaxies of later type show more pronounced 
differences between the undisturbed, blue model and the 
post-starburst phase of the burst model. Spectral transfor- 
mation into SO-type remnants is also slower for those galax- 
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ies since they overall have a relatively young stellar popula- 
tion which has to age first. 

Bursts with earlier onset times produce more luminous 
starburst and post-starbursts phases. The same holds true 
for increasing burst strength. The reason for both effects is 
our definition of the burst strength. The amount of stars 
formed during the burst is given by the available gas-mass 
at the onset of the burst multiplied with the burst strength. 
From there it follows naturally that earlier onset times cor- 
respond to a larger available gas-mass to be turned into stars 
as does a larger burst strength. 

Halting of star formation on the opposite leaves brighter 
galaxy remnants if it happened at a late stage, because until 
then star formation in the undisturbed galaxies leads to a 
continuous built-up of stellar mass. 

Bursts with longer decline time evolve slower in both 
colour and luminosity towards their postburst-phases and 
finally into SO-type galaxies. Also the luminosity increase at 
the peak of the burst decreases with longer decline times of 
the burst. 

In the case of a slow decline of SFR with t = 1 Gyr the 
evolution takes « 3.5 Gyr more time than in the case of a 
short decline with r = 0.1 Gyr. As mentioned above also the 
increase in luminosity is less prominent for the first case and 
colours do not reach as blue colours as in the latter case. 



3.3 Impact of bursts on the colour evolution 

During the burst and post-burst phases the U-B colour ex- 
periences the largest changes. Colours in the red part of 
the optical, e.g. V-I or R-J, also change slightly, while H-K 
colours hardly change at all. This can be explained by the 
stellar populations dominating each of those regions of the 
spectrum. The U and to a lesser degree the B to R-bands are 
dominated by young, hot stars which evolve rapidly with the 
declining star formation. NIR luminosities on the opposite 
are dominated by cooler and longer-lived stars and hence 
the increase in luminosity in these bands is roughly propor- 
tional to the increase in stellar mass, leaving the resulting 
colours essentially unchanged. 

The shape of the long wavelengths of the SED from I 
to K-band hence looks similar in all post-starburst galax- 
ies, while galaxies can have very different luminosities from 
-23.5 to -19 mag in the K-band. In general, all SEDs be- 
come fainter and redder after the burst with increasing age. 
For weak bursts galaxies can even become fainter than their 
undisturbed progenitor, since we assume that star formation 
declines exponentially to zero after the burst. 

In the case of SF truncation or termination the SED 
immediately becomes fainter than that of the undisturbed 
spiral galaxy, so there is no initial increase in luminosity 
characteristic for the galaxies with starburst. While the U 
to R part of the SED very quickly adopts the shape of an 
SO galaxy SED, the I to K part of the SED remains very 
similar to that of the undisturbed spiral galaxy. 



3.4 Progenitors and remnants 

In the upper part of Fig. O we show an Sa-type model 
with a burst of 50% at 6 Gyr and a decline time of 0.1 
Gyr with an age of 10 Gyr. We also show for comparison 
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Figure 3. SED of an Sa burst model (50% burst at 6 Gyr, decline 
time r = 0.1 Gyr) at 10 Gyr compared to the undisturbed E and 
SO models. 



the SED of an SO and elliptical model, both normalized 
to match the typical B -band luminosi t ies for their particu- 
lar type as taken from ISandage et al. 1 (|l985bl . Fig. 10) and 
ISandage et all (|l985al . Fig. 4). 

The SED of the 10 Gyr old Sa post-burst model is es- 
sentially indistinguishable from that of the SO galaxy. The 
luminosity difference between the SED for the E galaxy and 
the model SED is slightly larger than the 1 a luminosity 
ranges of SO a n d ellip tical galaxies of 1.5 mag as given by 
ISandage et"ai] (|l985bl . see also Table 0). The Sb and Sc 
model luminosities match the luminosities of observed SO 
galaxies within the 1 a range, but cannot reach the luminosi- 
ties of the typical ellipticals. The Sd model neither matches 
the SED of the E galaxy nor the SED of an SO galaxy in 
terms of luminosities. However, considering galaxy mergers 
and the involved increase in mass by a factor of « 2 for equal 
type mergers, all spiral model SEDs are able to match the 
overall SED luminosities of the SO galaxies approximately 4 
Gyr after their sta rbursts. This is in g ood agreement with 
results obtained bv lBicker et al.l (|2003 ). 

The lower part of Fig. [3] shows the magnitude differ- 
ence between the Sa burst model and the SO- and E-type 
model. We removed the offsets due to different stellar masses 
presented above by matching the I-band magnitudes of all 
models. This emphasizes the extremly similar shapes be- 
tween the post-burst Sa model and both E and SO. Colours 
longwards of R are essentially indistinguishable, only in the 
U-band are small differences. The E-type galaxy has even 
redder U-B colours than the Sa burst model in its postburst 
phase, while the SO-type model is slightly bluer in U-B. 

In Fig. |4]we subtracted the SED of the undisturbed Sa 
model from the SED of an Sa burst model with a 30% burst 
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Figure 4. Difference in magnitudes between tlie SEDs of an 
undisturbed and a burst Sa model at three different ages 11.5, 
12 and 13 Gyr. The undisturbed model was subtracted from the 
burst model (Sa, burst 30% at 11 Gyr, decline time 1.0 Gyr). 
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Figure 5. Difference in magnitudes between the SEDs of two 
burst models at three different ages 6.5, 7 and 8 Gyr. The Sd 
model (burst 50% at 6 Gyr, decline time 0.1 Gyr) was subtracted 
from the Sd model (burst 30% at 6 Gyr, decline time 0.1 Gyr). 



after 11 Gyr and a decline time of 1.0 Gyr at three different 
ages. This model shows one of the weakest bursts of the grid 
and does not even have strong Balmer absorption Unes. We 
only show the SEDs during and close to the E+A phase, i.e. 
at 11.5, 12 and 13 Gyr. Colours in the I through K-bands 
hardly change, only the luminosity changes towards brighter 
magnitudes as a result of the additional built-up of stellar 
mass during the burst. Towards shorter wavelengths the dif- 
ferences increase with both increasing age and decreasing 
wavelength, reaching up to 0.4 mag in the U-band 2 Gyr af- 
ter the burst. Since this model is one of the models with the 
weakest burst of our grid, the differences given for this Sa 
model are a lower limit. The differences for other models are 
larger and, for example, amount to up to 1 mag for the Sd 
model with the maximum peak value for H(5. We therefore 
conclude that post-starburst galaxies, even after very small 
bursts, can be distinguished from undisturbed galaxies on 
the basis of their SEDs. 

In Fig. [5] we compare SEDs of two different burst mod- 
els, both starting as Sd-type model and both encountering 
bursts with decline time 0.1 Gyr after 6 Gyr, but with dif- 
ferent burststrengths of 30% and 50%. 

We subtracted the SED of the weaker burst from the 
SED of the stronger burst, so that all values in Fig. [S] trans- 
late directly into magnitude differences between the two 
bursts in the sense that positive value mean that the stronger 
burst is brighter. The SEDs for ages of 6.5, 7, and 8 Gyr are 
shown. Although the stronger burst is brighter on average, 
there are only very small differences between their colours; 
identical colours in this representation would result in hori- 
zontal lines. 

At 6.5 Gyr the maximum difference in colours is 0.05 
mag and at 8 Gyr even less. Obviously some very similar 
models can not be clearly distinguished from each other by 
colours alone. However, the difference in luminosities is 0.4 
to almost 0.5 mag, which makes it possible to distinguish 
between models by luminosity. 

We therefore conclude that it is possible to distinguish 
E-j-A galaxies, post-starburst galaxies in a phase with still 
measurable [O ll] lines, and even post-starburst galaxies with 



weak H5-lines from undisturbed galaxies and, to a certain 
limit, from each other by carefully comparing their UV/U- 
opt-NIR SEDs. 

The advantage of this SED-based method is that only 
photometric observations have to be obtained while no spec- 
troscopic line measurements are necessary. Post-starburst 
galaxies with emission lines, which would be excluded by 
the classical E-j-A definition, can be detected and interpreted 
with this method, as well as post-starburst galaxies for which 
the burst was not strong enough to produce strong H5-lines 
with EW(H(5) ^ 5 A. It is also possible to distinguish be- 
tween most of the models, i.e. we can identify the most likely 
progenitor galaxy type and estimate the burst strength. We 
can thus investigate the full range of post-starbur st galax- 
ies without the selection effect we have identified in lPaper"^ . 
which excludes a significant number of post-starburst galax- 
ies from the E-l-A class. 

The U, B and V bands are the most important ones 
and hence should all be observed, while it is sufficient to 
h ave only o ne band in the near-infrared (NIR). The spectra 
in iPaper ll indicate that SEDs with additional filters in the 
UV, e.g. the FUV and/or NUV band from GALEX, would 
reveal even more significant differences between the SEDs of 
post-starburst galaxies. 



4 COMPARISON TO MORPHS DATA 

We will now use our models to study observed E-f A galaxies. 
We present the observational data and investigate them in 
terms of the Lick index H(5, luminosities, and colours. We 
also show how the best models can be found by comparing 
the SEDs of models and observations. 

We adopt a cosmological model with Ho — 
73kms~^ Mpc~\ Q m = 0.27, Qa = 0.73, and Zfo^m = 10 
(|Spergel et al.|[2007l ). The last parameter, the formation red- 
shift zform is the assumed redshift of galaxy formation Zform- 
We need this parameter to convert redshifts into galaxy ages 
at this redshift 

Tgal(z) = T(z) - T(zform) (1) 
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Table 3. Redshift as a function of galaxy age Tg^i for Ho = 
73kms-i Mpc^\ Qm = 0.27, Ha = 0.73, and zform = 10. 



T^gal [Gyr] Redshift z 
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where T(z) is the Hubble time, i.e. the age of the uni- 
verse, at redshift z. Tab. [3] shows the galaxy age Tgai as a 
function of redshift for the cosmology we chose. Galaxy ages 
and redshifts used in our models below are shown in bold 
font. 



4.1 The MORPHS E+A catalog 

For the subsequ ent analysis we use da t a from the MORPHS 
coUaboratioijj (|Poggianti et all 1 19991 : iDressler et"all 1 19991 : 
lOemlerlUgggf T They used the Wide Field Planetary Camera 
2 (WFPC2) onboard the Hubble Space Telescope (HST) to 
study the morphology of galaxi es in ten distant c lusters in 
a redshift range of z=0.37-0.56 (jSmail et al.ll 19971 '). 

The catalog lists redshift z, quality of the spectra, 
EW(H(5), luminosities, and colours for 88 possible E+A 
galaxies. The necessary spectroscopic data were obtained 
with the 200 inch (5.1 m) Palomar Hale Telescope (P200), 
the 4.2m William Herschel Telescope (WHT), and the 3.5m 
New Technology Telescope (NTT). Photometric data were 
taken with the 4-Shooter cameras on the 200 inch Palomar 
Hale Telescope and with HST using WFPC2. Ahhough the 
data have been taken with various telescopes, they were pro- 
cessed with the same data reduction methods to ensure ho- 
mogeneity. 

For most galaxies r, g— r and i in the Thuan & Gunn sys- 
tem (|Thuan fc Gunnlll976l . similar to the SDSS fihers, mag- 
nitudes are in AB system) of the P200 are given. Data from 
HST in I (F814W), R (F702W), B-I (F450W-F814W) and 
V-I (F555W-F814W) of the HST fiher system are given 
for only a few galaxies. To avoid uncertainties related to 
transformations between filter systems we therefore calcu- 
lated apparent magnitudes in all those filters by computing 
spectra first which are then convolved with the appropriate 
filter response curves. 



The data and catalogs can be found online at 
|http:/ / astro, dur.ac.uk/^ j irs/ morphs 



Figure 6. EW(H(5) versus redshift z of E-|-A galaxies observed 
by the MORPHS collaboration. 

4.2 Sample selection 

Galaxies with EW(H5) ^ 3 A and no significant emis- 
sion lines are classified as E+A galaxies. There is a sub- 
division into H(5-stronger E+A galaxies with H(5 ^ 8 A, 
called a+k galaxies, and H(5-weaker galaxies, named k+a, 
with 5A < EW(H(5) ^ 8A. We excluded aU galaxies with 
doubtful classification as w ell as spectra of qual ity 4 in the 
MORPHS grading scheme ([Dressier et al.|[l999l . where 1 is 
best and 5 is worst) to ensure reliable measurements of 
EW(H5). We also exclude galaxies with EW(H5) < 5 A, 
because H(5- values below this limit can easi ly be rea ched by 
non-poststarburst galaxies as we shown in IPaper j . Apply- 
ing these selection criteria leaves us with a final sample of 
44 E+A galaxies. 

4.3 Lick Index H5, Luminosities and Colours of 
Observed E+A Galaxies 

In Fig. [6] the EW(H(5) of all galaxies from the catalog are 
shown versus redshift. Filled circles indicate H5-strong a+k 
galaxies, open circles show H(5-weaker k+a galaxies. It can 
be seen that there is a group of galaxies at redshift z~0.4 
and another group at z ~ 0.55 corresponding to the redshifts 
of clusters studied by the MORPHS collaboration. The val- 
ues for EW(H5) range from 5 A (our selection criterion) to 
14.2 A. The measurements of E W(HJ) are reproducib le to 
±10% for H(5 > 5 A according to IDressler et al.1 |l99i). 

In Fig.[7l the EW(H(5) is plotted versus r- and i-band lu- 
minosities of the Thuan & Gunn system. Filled and open cir- 
cles again indicate a+k and k+a galaxies, respectively. The 
sys tematic errors for magn itudes and colours are estimated 
bv IDressler fc GunnI (|l992l ') to be of the order of 0.1 - 0.2 
mag. These observed values of EW(H(5) are to be compared 
with the evolution of our different model E W(H(5) during 
their post-burst phases as shown in IPaper j . Figs. 4-6. We 
recall that the maximum EW(H(5) reached by our models 
was around 8.5 A. Please note that the two galaxies with 
the highest values for EW(H5), located in the upper right 
part of the Fig. [71 are marked as questionable. We notice 
that there is an outlier in luminosity at r- and i-band lumi- 
nosities of 18.25 mag and 18.28 mag, respectively. However, 
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Figure 7. EW(H(5) versus Gunn & Tliuan r and i luminosities of 
E+A galaxies observed by the MORPHS collaboration. 



tliere is no note in tlie catalog tliat tliis galaxy has question- 
able values. If the luminosities and the E+A classification 
for this galaxy are correct, the brightness of this galaxy can 
only be explained by a strong previous burst in a merger of 
two previously very bright galaxies. 

Fig. |5] shows the HiS equivalent width EW(H5) as func- 
tion of the g-r colour in the Thuan & Gunn system. From 
the g-r colour plot we find that k+a galaxies are redder than 
a+k galaxies showing higher EW(H(5). This is in agreement 
with our models which show that galaxies with high EW(H5) 
can only have blue colours w hile galax ies with low EW(H5) 
can be either blue or red fsee lPaper J ). 

The provided sample with HST I and R luminosities 
as well as B-I and V-I colours is too small to make any 
statistical statements. The individual galaxies can still be 
analyzed, however, this is not the subject of this paper. 
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Figure 8. EW(H(5) versus colour for Gunn & Thuan g— r of E+A 
galaxies observed by the MORPHS collaboration. 



and convolved the spectra of our model grid with the Thuan 
& Gunn filter system and with the HST filter system. 

O ur GALEV models curr ently inc lude the H5 defini- 
tion of iTrager et al.l (|l998l . see iPaper j for details) which is 
slightly different from the definition of lBalogh et al.l (|l999l ) 
used by the MORPHS group. To ensure that this does not 
affect any of our results we measured EW(H5) of an A-star 
spectrum degraded to the resolution of the MORPHS spec- 
tra using both definitions. The result is that the maximum 
deviation amounts to 0.2 A, while the measurement unc er- 
tainties for the EW(H(5) is given bv lDressler et all (|l999l ) to 
be of the order of 10%, which corresponds to an uncertainty 
of 0.5 A for an EW(H5) of 5 A (the lowest value for EW(H5) 
of E+A galaxies). We therefore conclude that the EW(H5) 
as calculated by GALEV and the measured EW(H(5) of the 
MORPHS group are equivalent within the measurement un- 
cert ainties. 

iDressler et al.l (|l999l ) and lDressler et aL I (|2004 ) further- 
more compared the measurement of EW(H5) using the 
bandpass definition and direct line profile fits, the latter of 
which is used by the MORPHS group. They found that for 
weak but measurable cases, i.e. EW(H5) ^ 5A, the line- 
fitting and the bandpass methods are equivalent. For strong 
absorption lines, however, the bandpass definition tends to 
underestimate the EW, because in these cases the large pop- 
ulation of early A-type stars with relatively broad H5 that is 
necessary to reach these line strengths leads to a depression 
in the continuum bands, which in turns lowers the measured 
line strengths. This issue of measuring the li ne strength of 
H^ in an optimal way is discussed in length in IDressler et al.l 
(|2004 ). This explains why the observations in Figs.|6]|8]reach 
EW(H5) as high as 14 A, while the maximum value found 
in our model grid is only around 8.5 A. This does, however, 
in no way affect any of our results or conclusions. 



4.4 Methods of the Comparison between 
Observations and Models 

To compare our models with the observational data from 
the MORPHS catalog, we first have to ensure that EW(H(5) , 
colours and magnitudes are all given on the same system. 
For this purpose we implemented the Thuan & Gunn filter 
system of the Palomar Hale Telescope into the GALEV code 



4.5 Results from the Comparison between 
Observations and Models 

As we found in lPaperT . we can not reproduce the MORPHS 
observations by any of our models with SFR decline time 
r ^ 0.3 Gyr: Models with longer decline times either show 
remaining SF during the H(5 strong phases or are no longer 
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H(5 strong as soon as the SF has decreased sufficiently. How- 
ever, we will not exclude the models with a longer decline 
time to see how good the comparison of luminosities and 
colours between models and observations is when trying to 
find matching models. 

To compare MORPHS data at different redshifts with 
our models, the model spectra have to include cosmolog- 
ical and evolutionary effects as described in Sect. [2] Evo- 
lutionary (e-) corrections describe the fact that galaxies at 
redshift z — 0.4 are ~ 4 Gyr younger than local galax- 
ies at z = 0. Cosmological (k-) corrections on the other 
hand describe the bandpass shift due to cosmic expan- 
sion: The observed r-band magnitude (central wavelength: 
Ac ~ 6700 A) resembles approximately the rest-frame g- 
band (Ac = 5f00A fa 6700 A/(f 4-0.4)). To prevent fur- 
ther uncertainties from entering analysis we decided to not 
change the observed data, but rather include both e- and 
k-corrections into our models. 

At first we investigated which galaxy models are in their 
E-l-A phase in the redshift range between z = 0.36 — 0.56. 
Most models can already be excluded before any comparison 
with the data, because they are not in their H5-strong phase 
in the redshift range of interest. Of course, all galaxies with 
an onset of burst or truncation/termination at z ~ 0.35 can 
be excluded, because they are still undisturbed at 0.36 ^ 
z 0.56. 

The galaxy models with an onset at 2 ~ 0.93 can also 
be excluded because they have already passed through the 
H5-strong phase. ALL post-starburst galaxies with an onset 
at z ~ 0.44 reaching EW(H5) 5 A are possible. However, 
the galaxies with an onset at 2 = 0.44 only cover the low 
redshift part of the redshift range during their H^-strong 
phase. We can therefore already predict that galaxy models 
with an onset between z — 0.44 and z = 0.93 would fit the 
data much better. As examples a few galaxy models with 
an onset at 2: = 0.57 were calculated in addition to the grid 
explored so far. 

We furthermore chose to not only use colours but also 
absolute magnitudes in the subsequent analysis, hence not 
allowing galaxy masses as completely free parameter. The 
physical reason for this is that the properties during and 
after the starburst heavily depend not only on the amount 
of newly born stars, but also on the star formation history 
before the burst (see Sect. |3J. This spectral type depen- 
dence then translates into a mass-dependence, since galax- 
ies of each type only cover a narrow mass-range as de- 
rived from absolute m agnitudes in the local universe (e.g. 
ISandage et al.l [l985b| [all . This limitation hence rules out a 
number of otherwise possible scenarios, e.g. the formation 
of bright E-|-As by pure truncation of SF. 

In Fig. 191 141 we present luminosities and colours in 
various filters as function of redshift. We recall that our 
undisturbed model SFHs were chosen as to match after 13 
Gyr of undisturbed evolution the observed galaxies' average 
colours. Model masses were chosen as to yield the average 
B-b and lumi nosities observed for Sa, Sb, Sc and Sd galaxies 
(see iPaper j . Sect. 2.1). Observed undisturbed local galax- 
ies show 1 a ranges around those mean luminosities that 
need to be considered in the comparison between observa- 
tions and models. Very important to note is that on all the 
plots of EW(H(5), luminosities, and colours as a function of 
redshift the models evolve from the right-hand high redshift 
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Figure 9. HST F814W I-band luminosity versus redshift z. The 
symbols show the observed data. The box indicates the mean 
values with 1 cr standard deviations of the observation data. The 
Hues show some of the models with onset at 2: ~ 0.44 and ~ 0.93 
(top) and z ~ 0.57 (bottom). 



side towards the left-hand low red shift side in contrast to all 
the time evolution plots shown in lPaperT . 

Boxes indicate the mean values with 1 a ranges of the 
observational data in our sample. Shown are those models 
that have the maximum and minimum values of the grid and 
have their H5-strong phase in the given redshift range. We 
also show the undisturbed Sa and Sd models for comparison. 

At the top of Fig. [1 WFPC2 I-band luminosity is plot- 
ted versus redshift. The Sc galaxy with a burst of 70% and 
a decline time of 0.1 Gyr at ~ 0.44 has a luminosity, that 
matches quite well the two observed E-l-A galaxies at a red- 
shift of 2 ~ 0.4 (all models with onset 2 ~ 0.44 do so). 
Remember that the luminosities given in all plots are rep- 
resenting the average luminosity for each particular galaxy 
type, while true galaxies span a range in luminosities (see 
Tab. [21). Furthermore in the case of galaxies with bursts, 
those are induced by merging two galaxies, increasing the 
effective mass and hence luminosity by a factor of 2, cor- 
responding to a magnitude difference of 0.75 mag. Magni- 
tudes of E-l-A galaxies at z ~ 0.5 can be reproduced by 
earlier bursts. We show two examples of galaxies starting 
as Sc galaxies and encountering strong bursts of 70% with 
decline times of 0.1 and 0.3 Gyr. Both very well match the 
average i-band magnitudes observed at z « 0.5. 
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Figure 10. HST F555W V-band luminosity versus redshift z. 
The symbols show the observed data. The box indicates the mean 
values with 1 a standard deviations. The lines show some of the 
models with onset at z ~ 0.44. 



We also show for comparison the magnitude evolution 
of the undisturbed Sa and Sd models. The Sd is too faint by 
several magnitudes. The Sa-type model, however, is found 
to have approximately the right luminosity to matc h the ob- 
servations. However, as we have shown in IPaoe^ . Sa-type 
galaxies do not become H5 strong enough to be considered 
E-l-A galaxies, since they are lacking the gas reservoir to 
form a sufficient amount of new stars. Minor mergers of 
one early-type spiral and one late-type spiral could be a 
way to increase the gas-reservoir and turn early-type spirals 
into E+As. The detailed modelling of the range of possible 
unequal-type mergers, however, is beyond the scope of the 
present paper. 

In the bottom panel of Fig. |9] we compare bursts oc- 
curing in previously undisturbed Sd-type galaxies. Bursts 
start either at z = 0.57 or z — 0.44 and convert 70% of the 
then available gas mass into stars; decline times shown are 
0.1 and 0.3 Gyr. We manage to reproduce the luminosities 
of all E-l-A galaxies if we also account for the doubling of 
mass during the merger. Truncation scenarios, here shown 
for r = 0.1 Gyr and two different onset redshifts, however, 
are significantly too faint to match the observations. 

In Fig.[loiwe show the HST F555W (~ V) band mag- 
nitudes of all galaxies observed in this filter. Again we find 
that bursts of Sc galaxies with short decline times match the 
observations. Galaxies at higher redshifts can be explained 
by slightly earlier bursts with onset redshifts of 0.5 . . . 0.7. 
The only exception is one extremely bright galaxy with 
m_R ~ 18.2 mag or Mr = —23 mag using a distance modu- 
lus of m - M = 41.4 (Galaxy #18 of c luster CL0939-h4713 
in the catalog of Dressier et al.l (|l999l ) and # 224 in cata- 
log of lSmaU et all (|l997l )). This can only be explained by 
a massive merger of two above-average massive galaxies, 
or a complex merger of more than two galaxies, possibly 
a small, infalling group. Further evidence for this scenario 
is the complex shape (see Fig. Illfl showing significant sub- 
structure with multiple (tidal) tails as well as several nearby 
companion galaxies. However, on the basis of currently avail- 
able data we are not able to make more definite statements 
on the nature of this galaxy. 
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Figure 11. R-band image of CL0939+4713#18, taken with 
WFPC2. This galaxy shows a very irregular appearance with mul- 
tiple (tidal) tails sticking out of the main body. 
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Figure 12. HST F555W-F814W colour (~ V-I) versus redshift 
z. The symbols show the observed data. The box indicates the 
mean values with 1 a standard deviations. The lines show some 
of the models with onset at ^ ~ 0.57. 



Fig.[T2]shows the HST F555W-F814W (~ V-I) colour. 
We only show some examples for scenarios with burst begin- 
ning at redshift z ~ 0.57, because all observed galaxies lie at 
z ~ 0.55. The majority of the data points can be explained 
with these models. However, for some observed galaxies with 
red V-I colours, an even earlier onset at z ^ 0.6 has to be 
considered. 

In the top panel of Fig.[T3]the r-band luminosity in the 
Thuan & Gunn system is plotted versus redshift. All galaxies 
at z ~ 0.4 and z « 0.5 can be described by our models with 
strong bursts and decline times of 0.1 Gyr and 0.3 Gyr. The 
truncated Sd model is found to be too faint, in agreement 
with our findings for the HST filters above. 

The bottom panel of Fig. [13] shows the Thuan-Gunn 
i-band magnitude of all E-l-A galaxies. As found already for 
the g-band magnitudes we can reproduce the full spectrum 
of observed magnitudes. An important finding is that galax- 
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Figure 13. Thuan & Gunn luminosity in the r- and i-band ver- 
sus redshift z. The symbols show the observation data. The box 
indicates the mean values with 1 a standard deviations of the 
observation data. The lines show some of the models at z ~ 0.44 
and z ~ 0.57. 



ies which encounter a burst at redshift z « 0.57 are still 
bright enough to fit the data for those E+As at z « 0.4 
when we account for the 0.75 mag increase due to the dou- 
bling of mass in the merger. 

In both panels of Fig. [13] we also show the evolution of 
undisturbed Sa and Sd-type galaxies. Although the Sa-type 
model reproduces the required luminosities, both are not in 
an E+A phase at their particular age. Likewise in both plots 
we have one significant outlier at too bright magnitudes. 
This is the galaxy we have already discussed above. 

Summarising our discussion of model magnitudes as 
compared to observations we can conclude that many models 
can already be excluded on basis of magnitudes in individ- 
ual filters alone. Burst models in general reach the observed 
bright magnitudes, while models with SF truncation without 
previous starburst are too faint to match the observations. 

Fig.[l4]shows the g— r colour versus redshift. During and 
shortly after the burst our models are too blue to match the 
observations. However, models with strong bursts occuring 
at redshift z « 0.57 reach the red colours observed in E-|-A 
galaxies at redshift z « 0.44. This is a very important re- 
sult, since the time between those redshifts accounts to ~ 1 
Gyr (tgai(z = 0.57) = 7.3 Gyr and tgai(z = 0.44) = 8.2 
Gyr). After this time the SFR has declined sufficiently to not 
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Figure 14. Thuan & Gunn colour (g— r) versus redshift z. The 
symbols show the observation data. The box indicates the mean 
values with 1 a standard deviations of the observation data. The 
lines show some of the models at z ~ 0.44 and z ~ 0.57. 



show up in the [O II] lines anymore while EW(H5) reaches 
a maximum. We can hence fully reproduce all four factors 
determining E+A galaxies at those redshifts: 1) strong ab- 
sorption in H(5; 2) no emission in [O II]; 3) red colours 4) 
high luminosities. 

Galaxies observed in their E+A phase at redshift z ~ 
0.5 can be described by earlier bursts occuring at z « 0.75. 

Furthermore, only models with short decline times of 
^ 0.3 Gyr match the luminosities and colours, confirm- 
ing our results based on the EW([Oll]) described earlier. It 
seems that the constraints from luminosities and colours pro- 
vide very valuable complementary information to the studies 
of absorption lines and allow to discriminate between possi- 
ble and impossible scenarios. This proves that the compari- 
son of observed and model SEDs introduced in Sect. [3] is a 
good method to find adequate models for the observations. 



4.6 Comparison of Observed SEDs with Model 
SEDs 

In the following we will compare the full model SEDs to the 
observed SEDs. Since we have found earlier that our model 
magnitudes match the average observed magnitudes for the 
E+A galaxies, we normalize all SEDs to a r-band magni- 
tude of 0, so that our SEDs essentially consist of the g-r and 
i-r colours. We do not show the HST magnitudes, since the 
majority of galaxies in the MORPHS catalog were only mea- 
sured in the Thuan-Gunn system. Fig.[T5]shows the observed 
colour range of all a+k and k+a galaxies as grey-shaded ar- 
eas. In the top panel of this Fig. we also show three dif- 
ferent model-SEDs based on Sc-type progenitor galaxies for 
redshift z = 0.4. In accordance with our analysis of the indi- 
vidual colours we find that the model with burst at z = 0.44 
is still too blue to explain the observations. Models with ear- 
lier bursts ai z — 0.57 lie within the observed colour range, 
with the shorter burst-duration of 0.1 Gyr providing a bet- 
ter match than the model with the slightly longer duration 
of 0.3 Gyr. 

The lower panel of Fig. [15] compares the observed data 
with Sd-progenitor models. Both models with burst redshift 
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Figure 15. Comparison of full SEDs. The range spaned by the 
observed a+k (k+a) galaxies are shown as dark (light) shaded 
grey region, In the top panel we show SEDs for three burst models 
with Sc progenitors, in the lower panel we show four models with 
Sd-typc progenitor galaxies. 

of z = 0.57 and burst durations 5C 0.3 Gyr lie within the 
observed range. The model with the longer burst duration 
of 1 Gyr is still too blue at redshift z = 0.4 for which the 
SEDs are shown. For comparison we also show an SED of an 
earlier burst at z — 0.9. Although this matches the observed 
SED quite well, the galaxy model is no longer in an E+A 
phase at the observed redshift z = 0.4, so that this model 
can be excluded as well. 



5 WHERE DOES THE REMAINING GAS GO? 

One question remains yet untouched: The fate of the re- 
maining gas that was not converted into stars during the 
burst. In the following we will show that the energy input 
by supernovae is sufficient to remove this gas. Integration 
of an Salpeter IMF with mass- limits of 0.1 and 100 Mq 
yields that ~ 10% of the gas-mass is converted into stars 
with M > 8M0, that end their life as core-collapse super- 
novae. We will base our example on a burst-model starting 
as Sc-type model with total mass of 1.4 x lO^'' Mq and en- 
countering a strong (70 %) burst at an age of 7.5 Gyr. This 
burst forms a total stellar mass of 4.8 x 10^ Mq. We further 
assume a typical mass of IOM0 for stars with M > 8 Mq 



and that each resulting supernova produces a total energy of 
10^^ erg (|Burrows et al.lll995l : lHamuvll2003l '). All supernovae 
formed during the burst hence produce a energy output of 
« 5 X 10^® erg. To derive the potential energy needed to 
remove the remaining 2 x 10^ Mq gas we use the total mass 
of the galaxy and assume a typical radius of 3 kpc. This re- 
sult in a binding energy of ~ 7.3 x 10^^ erg, or 14% of the 
energy produced by the supernovae. This fraction is in good 
agreement wit h earlier results o f 0.09 l|Thornton et al.lll998l ) 
and 0.1 - 0.2 (|CoIe et al.lll994h . 



6 SUMMARY AND OUTLOOK 

Following the first part of our study presented in IPaper j 
we present Spectral Energy Distributions for yet undis- 
turbed, pre-starburst galaxies, galaxies with ongoind star- 
bursts and post-starburst galaxies. We showed that com- 
paring the SEDs of post-starburst galaxies to those of undis- 
turbed spiral, SO and E galaxies as well as among each other 
is a sufficient method to select and investigate post-starburst 
galaxies. The interplay between the burst strength, the de- 
cline time, the onset of the burst/halting of SF and the pro- 
genitor galaxy type shapes the SED in the following way: 

(i) The later the progenitor galaxy type, the more obvi- 
ous is the difference between the SED in the post-starburst 
phase to the SED of the corresponding undisturbed galaxy 
model. Late- type spirals furthermore evolve more slowly to- 
wards SO spectral types than early-type spirals. 

(ii) The earlier the onset of the burst the more luminous 
is the model SED in its starburst and post-starburst phases. 

(iii) The higher the burst strength, the more luminous is 
the model SED throughout its starburst and post-starburst 
phase. 

(iv) The longer the decline time, the longer takes the time 
evolution of luminosities and colours of the SED to evolve 
first into the SED of a post-starburst galaxy and later into 
the SED of an SO galaxy. 

We have shown that we are able to distinguish E-l-A 
galaxies, post-starburst galaxies in an early phase after 
a burst with still measurable [0 11] lines, and even post- 
starburst galaxies with weak H(5-lines from undisturbed 
galaxies and, to a certain limit, from each other by care- 
fully comparing their SEDs. 

The advantage of this method is that only photomet- 
ric observations have to be made and no line measurements 
are necessary. Post-starburst galaxies with emission lines, 
which are excluded by the classical E-f A definition can be 
detected and interpreted, as well as post-starburst galaxies 
for which the burst was not strong enough to produce strong 
H5-lines with EW(H(5) > 5 A. It is possible to distinguish 
between most of the models i.e., we can identify the progeni- 
tor galaxy type and estimate the burst strength. This means 
we can investigate the full range of post-starburst galaxies 
without the selection effect which excludes a number of post- 
starburst galaxies from the E-f A class (see lPaper J) . 

The most important bands crucial for this SED anal- 
ysis are U, B and V, while it is suffic ient to h ave only one 
band in the NIR. Our model spectra in lPaperT indicate that 
SEDs with additional filters in the UV, e.g. from GALEX, 
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would reveal even more subtle differences among the SEDs 
of different types of post-starburst galaxies. 

In the second part of this paper we compared our model 
results to the MORPHS catalog containing 88 E+A galaxies 
in the redshift range z=0.36— 0.56. We constructed a grid of 
models spanning a wide range of progenitor galaxies, burst 
redshifts, burst decline times and burst strengths. All our 
models fully include both cosmological and evolutionary cor- 
rections as function of redshift and galaxy type. 

We find that k-l-a galaxies are redder than a-|-k galaxies 
with higher EW(H5). This is in good agreement with our 
models which show that galaxies with high EW(H5) can 
only have blue colours while galaxies with low EW(H5) can 
either be still blue or already red. 

From our comparison of model and observational data, 
we find that only galaxies with strong bursts ~ 70% 
and relatively short decline times ^ 300 Myr can ex- 
plain the observed colours and magnitudes. This is in good 
agreement with earlier studies (iPoggianti fc Barbard 1 19961 : 
iBalogh et all Il997l . l2005l : lYan et al.ll2008h which also found 
indications for recent bursts forming a significant fraction 
of the observed stellar masses. Weaker bursts do not reach 
the required H(5 strengths, while models with SF truncation 
instead of a starburst are significantly too faint, con firm- 
ing earlier findings from IChristlein fc Zabludofj (|2004l ) who 
found that simple fading without an increase of the bulge 
luminosity does not reproduce observations of SO galaxies. 

We hence conclude that a large fraction of the parame- 
ter space can easily be excluded based on SEDs alone. This 
is remarkable since our SEDs only consist of three filters and 
only cover a short wavelength range of the spectrum. 

We could also show that supernovae being formed as a 
consequence of the burst can produce enough energy to re- 
move the remaining gas that was not turned into stars from 
the galaxy, explaining why most observed E-|-A galaxies do 
not contain significant amounts of cool gas. 

The predictions from our models can be used to col- 
lect more and better data for the SEDs of post-starburst 
galaxies. With better sampled SEDs, in particular with UV 
filters and one NIR band included, and a su itable SED fit- 
ting tool (e.g. AnalySED (|Anders et al.ll2004l 'l or GAZELLE 
(KotuUa & Fritze, in preparation)) it will become possible 
to improve the method outline here and to better constrain 
not only progenitor galaxies and transformation scenarios, 
but also limit the range of possible burst ages and burst 
fractions. 

Investigating SEDs of all post-starburst galaxies 
promises to lead to a better understanding of the transfor- 
mation processes of galaxies in high density environments. 
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